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1.  INTRODUCTION 


Much  effort  in  recent  years  has  been  devoted  to  formulating  an 
accurate,  yet  tractable,  theoretical  procedure  for  analyzing  the  optical 
spectra  of  triply  ionized  lanthanides  in  crystals.  Owing  to  the  many 
levels  in  the  f^  ground  configurations  of  these  ions  and  to  the 
relatively  few  low-lying  levels  whose  Stark  splittings  are  observed,  it 
is  undesirable  and  often  infeasible  to  require  simultaneous 
diagonalization  of  the  free-ion  Hamiltonian,  H,  and  the  crystal  field 
Hamiltonian,  Hjj.  It  is  preferable  first  to  diagonalize  H,  the  larger 
interaction,  and  then  to  calculate  the  crystal  field  splittings  by 
diagonalizing  Hy  in  the  low-lying  levels.  This  report  examines  this 
procedure,  specifically  for  calculating  Pr^'*’,  Tm^'^,  Nd3  + , and  Er3  + . 

Theory  cannot  yet  determine  the  free-ion  Hamiltonian  ab  initio,  but 
its  form,  H(P.),  is  adequately  defined  by  a set  of  P.  that  determine  the 
strengths  of  the  Coulomb,  the  spin-orbit,  and  other  interactions. 
Values  for  these  parameters  are  defined  as  minimizing  the  difference 
between  theoretical  free-ion  centroids  evaluated  in  a basis  of  inter- 
mediate coupled  eigenfunctions  and  experimental  free-ion  centroids. 
Since  few  experimental  data  are  available  for  the  f^'  configurations  of 
free,  triply  ionized  lanthanides,^*^  the  P.,  in  practice,  are  determined 
to  give  a best  fit  to  Stark-level  centroids  in  aqueous  solution  or  in 
crystals.  Unfortunately,  the  crystal  field  affects  centroid  energies 
so  that,  even  if  all  Stark  components  of  a given  J-multiplet  are 
determined  experimentally,  the  weighted  centroid  is  not  the  free-ion 
energy  resulting  if  the  crystal  field  were  "turned  off."  As  a second 
problem,  even  the  most  elaborate  free-ion  Hamiltonian  often  cannot  give 
an  adequate  fit  to  the  J-multiplet  centroids  once  they  are  determined. 
The  seven-parameter  Hamiltonian  defined  by  Rajnak  and  Wybourne^'^  and 
used  extensively  by  Carnall,  Fields,  and  Rajnak  (CFR),**  for  example, 
often  gives  centroid  discrepancies  of  100  cm*^. 

This  report  describes  a method  that  circumvents  the  first  of  these 
difficulties.  The  method  is  based  on  the  fact  that,  in  perturbation 
theory,  changes  in  wave  functions  give  second -order  changes  in  energy. 
Conversely,  known  small  changes  in  energy  determine  the  perturbing 
Hamiltonian  to  first  order.  Thus,  any  reasonable  set  of  P^  presumably 


^J.  Sugar,  Phys.  Rev.  Lett,,  (1965),  731. 

^N,  Spector  and  J.  Sugar , J,  Opt.  Soc . Am.,  6^  (1976),  436. 

^V.  Kaufman  and  J.  Sugar,  J.  Opt.  Soc.  Am.,  6^  (1976)  439. 

‘*W.  T.  Carnall,  P.  R.  Fields,  and  K.  Rajnak,  J.  Chem.  Phys.,  49 
(1968),  4412-55. 

^K,  Rajnak  and  B.  G.  Wybourne,  Phys.  Rev.,  1 32  (1963),  280. 

^B.  G.  Wybourne , Spectroscopic  Properties  of  Rare  Earths, 

Interscience  Publishers,  New  York  (1965). 


determines  an  adequate  set  of  J-multiplet  wave  functions  and  reduced 
matrix  elements,  even  though  the  energies  they  determine,  E^,  may  be 
incorrect.  An  appropriate  subset  of  these  "free-ion"  J-multiplet  wave 
functions  (usually  10  to  15  of  the  lowest  lying)  may  then  be  chosen  as  a 
basis  to  calculate  Stark  splittings,  and  the  energies,  E^ , of  these 
multiplets  may  be  treated  as  parameters  (along  with  the  crystal  field 
parameters,  in  Hj{)  to  be  varied,  as  suggested  by  Margolis,^  for 

obtaining  a fit  between  theoretical  and  experimental  Stark 

levels.®'^  The  resultant  E*,  which,  in  general,  differ  (by  AEj)  from 
the  initial  centroids,  may  finally  be  used  to  determine  new 
+ AP , , as  suggested  by  Wong,^®  by  minimizing  the  function 

Q = I .ri.  (3E./3P, ) AP.  - AE.]2  (1) 

with  respect  to  AP^.  If  the  APj^  determined  in  this  way  are  not  large, 

then  the  linearity  assumption  in  Wong's  expression  is  valid,  and  also  it 

is  unnecesarry  to  rediagonalize  the  free-ion  Hamiltonian  to  determine 

new  reduced  matrix  elements.  These  results  suggest  that  it  is  desirable 

to  perform  the  initial  diagonalization  of  H(P^)  with  P,  that  pertain  to 

some  convenient  host  and  not  to  the  true  free  ion  since  variations  in 

the  P.  from  host  to  host  are  expected  to  be  less  than  the  change  in  the 

P.  from  the  free  state  to  some  host.'^^ 

1 

Accordingly,  we  have  calculated  theoretical  energies,  derivatives, 
and  reduced  matrix  elements  for  the  f^,  f*^,  f^,  and  f^^  conf igurations , 
respectively,  of  Pr3+,  Tm3+,  Nd^'*',  and  Er3'*‘  using  the  aqueous  solution 
free-ion  parameters  of  CFR.**  Wong^®  determined  the  3E^/3P^  for  Nd^"*”  and 
Er^"'"  for  P.  = {E^,  E^,  E^,  by  a finite  differencing  method  involving 
five  separate  diagonalizations  for  each  ion.  We  have  extended  and 
improved  on  his  work  by  obtaining  derivatives  also  with  respect  to 


‘*W.  T.  Carnall,  P.  R.  Fields,  and  K.  Rajnak,  J.  Chem.  Phys.,  ££ 
(1968) , 4412-55. 

S.  Margolis,  J.  Chem.  Phys.,  £5  (1961)  1370. 

Karayianis,  D.  E.  Wortman,  and  H.  P.  Jenssen,  J.  Phys,  Chem. 
Solids,  J7  (1976),  676. 

^P.  Grunberg,  S.  Hufner,  E.  Orlich,  and  J.  Schmidt,  Phys.  Rev.,  184 
(1969),  290. 

Y.  Wong,  J.  Chem.  Phys.,  3^  (1961),  544. 

Karayianis  and  C.  A.  Morrison,  Rare  Earth  Ion-Host  Crystal 
Interactions.  2.  Local  Distortion  and  Other  Effects  in  Reconciling 

Lattice  Sums  and  Phenomenological  B , Harry  Diamond  Laboratories 
TR-1682  (January  1975),  table  II, 


Rajnak  and  Wybourne's^  higher-order  interactions  of  strengths  a,  B/ 
and  Y and  by  using  Feynman's  theorem'^  to  obtain  derivatives.  Feynman's 
theorem,  expressed  by 


3E  ./3P.  = <11/ . 3H/3P,  U .> 

: 1 i'  3 


(2) 


where  represents  the  jth  eigenfunction  of  H,  has  the  advantages  of 
having  "^greater  inherent  accuracy  and  of  requiring  only  one 
diagonalization  per  ion.  To  perform  the  calculations,  electrostatic 
matrices  corresponding  to  the  Russell-Saunders  basis  states  given  by 
Nielson  and  Roster  were  hand  keypunched  from  their  book  of 
tables.  The  v and  U reduced  matrix  elements  that  relate, 
respectively,  to  spin-orbit  and  spherical  harmonic  reduced  matrix 
elements  were  calculated  from  a computer  tape  containing  coefficients  of 
fractional  parentage  kindly  provided  by  H.  M.  Crosswhite  (Argonne 
National  Laboratory) . By  choosing  to  diagonalize  H with  the  aqueous 
solution  parameters  of  CFR,**  we  were  able  to  check  our  computer  software 
procedures  by  comparing  our  centroid  energies  with  those  of  CFR. 

A summary  of  the  computational  procedures  and  formulas  is  given  in 
section  2.  Section  3 gives  tables  of  the  centroid  energies  and 
derivatives  and  notes  some  discrepancies  between  our  Nd^'*'  and 
Er^'*'  energies  and  those  reported  by  CFR.*^  Section  4 gives  tables  of 
reduced  matrix  elements  of  the  spherical  harmonics  and  discusses  their 
use  in  crystal  field  calculations.  Also  in  section  4,  some  measure  is 
given  of  the  sensitivity  of  theoretical  Stark  levels  to  the  free-ion 
parameters  and  wave  functions  for  Nd^”^  and  Er3+  in  CaW04 . 


2.  CALCULATIONS 

To  obtain  free-ion  wave  functions,  matrix  elements  of  the 
Hamiltonian 


H = He  + Hg^  + ol2  + 6G2  + YR7 


(3) 


‘*W.  T.  Carnall,  P.  R.  Fields,  and  K.  Rajnak,  J.  Chew.  Phys.,  49^ 
(1968) , 4412-55. 

^K.  Rajnak  and  B.  G.  Wybourne , Phys.  Rev.,  132  (1963),  280. 

Merzbacher , Quantum  Mechanics , John  Wiley  and  Sons,  Inc.,  New 
York  (1961),  396. 

W.  Nielson  and  G.  F.  Foster,  Spectroscopic  Coefficients  for  the 
p^ , d'^ , and  f^'  Configurations,  The  MIT  Press,  Cambridge,  MA  (1963). 
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) 


(4) 


were  taken  in  a given  J subspace  spanned  by  the  basis  functions* 


JMLSo)> 


I <L(m)S(M  - m)  J (M' > Lm  S M-m  cj> 
m ' 


where  < jj  (ini ) j2  (m2)  | j (m)  > is  a Clebsch-Gordon  coefficient,  and  the 
I LSw>  are  the  Russell-Saunders  states  defined  in  Nielson  and 
Koster.^^  The  additional  quantum  number  oj  represents  the  pairs  and 
triplets  of  quantum  numbers  ujU2  and  WjW2W3  given  for  each  state  in 
Nielson  and  Roster  and  corresponding,  respectively,  to  the  G2  and 
R7  operators,^ 

Since  H is  a spherical  tensor  of  rank  zero  in  J space,  matrix 
elements  of  H in  a given  J subspace  are  independent  of  M.  If  the 
quantum  number  M is  suppressed,  the  electrostatic  matrix  elements  are 


1 


<JL'S'u'|h  |JLSa)>  = 6 ,6  ,<LSu)'|h  |lSw>  ; (5) 

£ LL  SS  E 

the  spin-orbit  matrix  elements  for  N equivalent  electronics  of  orbital 
angular  momentum  i are 

<JL'S'w'  lH^^jjLSio>  = -id  a + 1)  [£])*’ 

X WCIS'LJ;  SL')  (L'S'w*  I I |lS(jj)  ; (6) 

and  the  matrix  elements  of  the  last  three  operators  in  equation  (3)  are 


<JL'S'u'  h.  JLSw>  = 6,, ,5^^. 5 <h . > 
'1'  LL'  SS'  oju'  1 


(7) 


where 


= aL(L  + 1);  8G2(UiU2);  YH7(wiW2W3) 


(8) 


*Note  the  order  of  L and  S in  this  definition  since  it  does  make  a 
difference  in  the  matrix  elements . 

G,  Wybcurne , Spectroscopic  Properties  of  Rare  Earths, 
Interscience  Publishers,  New  York  (1965) . 

^ W.  Nielson  and  G.  F,  Foster,  Spectroscopic  Coefficients  for  the 

p , d^ , and  f^  Configurations,  The  MIT  Press,  Cambridge,  MA  (1963). 

^‘^M,  Rotenberg , R.  Bivins,  N.  Metropolis , and  J.  K.  Wooten,  The  3-j 
and  6-j  Symbols,  Technology  Press,  Cambridge,  MA  (1959) . 

M.  Rotenberg,  R.  Bivins,  N.  Metropolis , and  J.  K.  Wooten,  The  3-j 
and  6-j  Symbols,  Technology  Press,  Cambridge,  MA  (1959).  In  terms  of 
their  tabulated  3-j  symbols,  <j\(mi)j2(m2)\j(m)>[j]  = (-)~jdJ2~^ 

(h  72  j) 
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I 
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The  symbol  W(abcd;ef)  in  equation  (6)  is  a Racah  coefficient^  and 
[x]  = 2x  + 1.  Values  for  the  electrostatic  matrix  elements  and  the 
V n 1)  reduced  matrix  elements  are  given  in  Nielson  and 

Koster;^^  however,  we  keypunched  only  the  former.  A generalized 
defined  by 


(L'S'u)'||  I |lSw)  = N{s(s  + 1)  [s]  [L]  [L-]  [S]  [S']  } 


X E—  W(L'Lki;X.L)W(S'SKS;sS) 


X E-  F(L'S'w’N;LS(L  N-l)F(LSwN;LSiL  N-1)  , (9) 

u 

where  s = 1/2,  allows  one  to  calculate  a variety  of  reduced  matrix 
elements  once  the  coefficients  of  fractional  parentage, F,  are 
available.  From  H.  M.  Crosswhite's  (Argonne  National  Laboratory) 
computer  tape  containing  the  F for  the  configurations,  we  calculated 
the  V , as  well  as  the  reduced  matrix  elements  of  the  spherical 

harmonics,  C^,  according  to 

I I - ‘ss'l-TTiTInil}'’ 


X <i  (O)k(O)  1 (0)  >W(JkSL'  ; J'L)  (L'So)'  I |V 


LSw)  , * (10) 


where  H = 3. 

These  reduced  matrix  elements  are  necessary  for  the  crystal  field 
calculations  discussed  in  section  4.  The  reduced  matrix  elements,  in 
general,  are  defined  by 


<j 'm' jm>  = < j (m)  k (q)  I j ' (m' ) > [ j ’ ] (j’(|Tj^|(j) 


(Ok)  (k) 

*This  V is  related  to  the  U in  Nielson  and  Koster  by 

(\\v(0k)\\)  = (scs  + 1)  (S]]^(  \ u^-‘>\  \ ) . 

W.  Nielson  and  G.  F.  Koster,  Spectroscopic  Coefficients  for  the 
p'<,  d'^ , and  f Configurations , The  MIT  Press,  Cambridge,  MA  (1963). 

Rotenberg,  R.  Bivins,  N.  Metropolis , and  J.  K.  Wooten,  The  3-j 
and  6-j  Symbols,  Technology  Press,  Cambridge,  MA  (1959).  In  terms  of 

a+b+c+d (a  b e \ 

their  tabulated  6-j  symbols,  W(abcd;ef)  = (-)  \d  c f ) ' 

^^B.  R.  Judd,  Operator  Techniques  in  Atomic  Spectroscopy , McGraw-Hill 
Boole  Co.,  Inc.,  New  York  (1963),  166-192. 
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Values  for  the  eigenvalues  of  the  G an>-.  R7  operators,  GaCuju?)  and 
R7(wiW2W3),  are  given  in  table  I,  which  is  a copy  of  Wybourne's^  tables 
2-6  and  2-7. 

Taking  matrix  elements  according'  to  equations  (5)  , (6)  , and  (8) 

between  all  | LSo)>  states  that  satisfy  |l-s|^  J £ |l  + s|,  one  may  then 

choose  values  for  E^,  (in  ^ 

diagonalize  the  subspace  to  obtain  free-ion  wave  functions, 

1j(LSu))>  = S— -A— ,(LSw)  |jLSw>  (12) 

LSo)  LSu 

The  orbital  and  spin  angular  momenta,  L and  S,  in  |j(LScj)>  are  not  good 
quantum  numbers  because  of  mixing  by  the  spin-orbit  interaction, 
and  w is  not  a good  quantuln  number  because  of  mixing  by  both  the 
electrostatic  and  spin-orbit  interactions.  Thus,  the  "almost  good 
quawtum  numbers"  (LSoj)  in  parentheses  represent  values  in  the  limit  of 

the  off-diagonal  parts  of  and  H50  approaching  zero.  There  is  surh 

large  mixing  in  many  cases,  however,  that  the  assignment  of  LSoj  valuta 
to  a given  free-ion  wave  function  is  arbitrary.  The  method  that  we  have 
choosen  to  make  assignments  is  described  in  section  3. 


TABLE  I.  VALUES  OF  G2(uiU2)  AND  R;(wiW2W3) 


U1U2 

I2G2 

w^W2W3 

5R7 

00 

0 

000 

0 

10 

6 

100 

3 

1 1 

12 

110 

5 

20 

lA 

111 

6 

2' 

21 

200 

7 

22 

30 

210 

9 

30 

2A 

211 

10 

3' 

32 

220 

12 

AO 

36 

221 

13 

222 

15 

G.  Wybourne,  Spectroscopic  Properties  of  Rare  Earths, 
Interscience  Publishers,  New  York  (1965). 


3.  FREE-ION  ENERGIES  AND  DERIVATIVES 


Tables  II  to  V give  the  assignments,  energies,  and  derivatives  for 
each  of  the  free-ion  multiplets  of  Pr^"',  Tm3+ , Nd3+,  and  Er^"^, 
respectively.  The  arbitrariness  in  assigning  a particular  term  to  some 
highly  mixed  multiplet  was  resolved  by  scanning  the  matrix  of 
eigenvectors  for  the  largest  magnitude  coefficient,  assigning  the  row 
designation  of  that  coefficient  to  the  column  eigenvector,  eliminating 
that  row  and  column  from  further  consideration,  and  repeating  the 
procedure  until  the  final  eigenvector  was  assigned  by  default.  The 
square  of  the  designating  coefficient  is  printed  under  the  "PCT.  PURE" 
column  of  these  tables  and  indicates  the  percentage  of  the 
Russell-Saunders  designation  present  in  the  free-ion  eigenfunction. 
This  procedure  is  not  entirely  satisfactory  (although  it  results  in  a 
unique  one-to-one  assignment  for  a given  set  of  free-ion  parameters) 
because  a situation  may  arise  where  the  Russell-Saunders  designation  is 
not  the  largest  component  present  in  an  eigenvector.  The  possible 
occurrence  of  this  situation  is  for  eigenvectors  with  a PCT.  PURE 
assignment  of  less  than  50  petcent.  The  ^Ggy2(D  level  of 
Er^'^  (table  V)  is  an  example  of  this  occurrence.  It  contains 
24.23  percent  of  ^Fg/2»  but,  since  that  designation  had  already  been 
assigned  to  the  level  at  15,136  cm~^  which  contains  57.64  percent 
of  next  largest  percentage,  19.06  percent  of  '-09^2(1)' 

assigned. 

The  derivatives  with  respect  to  a yield  an  added  dividend 
since  3E/3a  = <L^>  according  to  equations  (2)  and  (3) . Thus,  the 

effective  orbital  angular  momentum  of  an  eigenvector  given  by 

= -1/2  + (1/4  + 3E/3a)*^  (13) 

gives  an  indication  of  the  L-mixing  induced  by  the  spin-orbit 

interaction.  For  the  ^Ggy2(l)  level  of  Er^'*',  L = 4.32,  indicating 

that  good  percentages  of  Hgy2  and  19/2  components  also  are  i^resent. 

Table  VI  gives  the  free-ion  parameters  for  these  ions  in  aqueous 
solution  determined  by  CFR.**  Our  calculated  energy  levels  agree  with 
theirs  for  Pr^"^  and  Tm^*.  For  the  ^Hjjy2(2)  level  of  Nd^*,  we  get 
16,043  (adding  130  cm‘3  to  our  value  in  table  IV  to  adjust  for  their 
ground  state)  compared  with  their*^  listed  16,026  cm"^.  In  Er^"^,  we  get 

47,811  and  36,572  cm"^  for  the  ^Hgy2(l)  and  ^Hgy2(2)  levels, 

respectively  (adding  109  cm'^  to  our  values  in  table  V) , while  they" 
list  47,822  and  36,566  cm"^ — differences  of  11  and  6 cm"l,  respectively. 
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T.  Carnall,  P,  R.  Fields,  and  K.  Hajnak, 
(1968),  4412-55. 


J.  Chew,  Phys.,  4° 
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TABLE  VI.  FREE-ION  PARAMETERS  FOR  TRIPLY  IONIZED  Pr,  Nd , Er, 
AND  Tm 


C ■'»  i*' 

IUO.15  21.255  -799-9't 

1312.9 

881. 58  0.5611  -117-15 

1321.3 

877.25  0.336  -112.52 

1211.3 

906.00  0.0  0.0 

0.0 

2380.7  18.317  -509.28 

619.71 

2399.7  23. 137  -626.01 

883-01 

2628.7  11.677  -631.79 

0.0 

‘‘for  aqueous  solution.  Sec  i T.  Carnall  et  al,  J.  Chcm.  I’hys . , 
^ (19b8),  ■4‘JlJ-55. 

^for  CaWOi.. 

^For  LaCl].  See  F.  Y.  iVonq , J.  Chan.  Fhys . , 35  il9ul),  5qq. 


Eight  other  levels  differ  in  energy  from  2 to  4 cm“^  in  Er^'*'.  We  do  not 
know  the  source  of  these  discrepancies,  and  we  solicit  results  from  an 
independent  calculation.  One  further  contact  was  made  with  Wong's 


Nd^'*'  calculations^^  where  we  found  agreement  to  0.6 


except  for 


his  level  at  19,288.93  cm''^.  We  calculate  this  energy  as 
19,228.9  cm"^,  strongly  indicating  a misprint  in  Wong's  table  that  would 
be  easy  to  miss  in  proofreading.  We  also  agree  with 
his  9E/3^  calculations  for  all  his  levels  with  spot  checks  with 
his  3E/3Fj^  values,  converting  them  to  our  3E/3e'^. 


4.  Wd  AND  Er  IN  CaWO,, 


Using  the  procedure  described  in  section  1,  we  determined'”  free-ion 
centroids  for  Nd^"'  and  Erl+  and  crystal  field  parameters  in  CaW04  by 
fitting,  respectively,  to  51  (of  64)  and  38  (of  48)  identified 
experimental  Stark  levels.  The  Nd^'*'  centroids  for  the  lowest  14 
multiplets*  were  fit  by  using  the  derivatives  of  table  IV  to  determine 

*Actually,  the  multiplet  is  lower  than  the  multiplct 

that  was  the  highest  multiplet  diagonalized  (see  table  IV) . However, 
the  number  of  Stark  levels  in  ^k'13/2  exceeded  the  available  core 
capacity  in  the  IBM  7094  computer  used  to  perform  these  calculations,  so 
the  '*Gij2  ttmltiplet  was  used  instead. 

Y.  Wong,  J.  Chem.  Phys.,  3^  (1961),  544. 

^^Donald  E.  Wortman , Clyde  A.  Morrison,  and  Nick  Karayianis , Rare 
Earth  Ion-Host  Lattice  Interactions.  5.  Lanthanides  in  CaWO^ , Harry 
Diamond  Laboratories  TR-1794  (June  1977). 


the  free-ion  parameters  for  CaW04  given  in  table  VI,  line  3.  Table  VII 
gives  theoretical  centroids  by  using  aqueous  parameters,  theoretical 
centroids  by  using  CaW04  parameters,  and  experimetal  centroids  for 
comparison.  The  listed  centroids  corresponding  to  CaW04  parameters  were 
obtained  by  Wong's  procedure  by  using  our  derivatives,  yet  they  differ 
only  by  0.27  cm"^  rms  from  the  exact  values  obtained  by  rediagonalizing 
the  free-ion  matrices.  The  new  eigenfunctions  and  reduced  matrix 
elements  obtained  by  rediagonalizing  were  then  used  to  recalculate  the 

TABLE  VII.  Nd3+  MULTIPLE!  CENTROIDS  (cm'O 


Mul t ip)et 

Carnal ) 

^ a 

parameters 

Best  fit^ 
parameters 

.c 

Experimental 
(CaW04  centroids) 

"1 

9/2 

8L 

219 

222 

4| 

11/2  1 

1961 

2081 

2089 

“1 

13/2 

3958 

4063 

4066 

15/2 

goig 

6120 

6117 

4p 

3/2 

11^81  ' 

1 

11488 

11453 

"F 

5/2 

12528 

12528 

12503 

(2) 

9/2 

12692 

12632  1 

1 

12597 

"S 

3/2 

13'tl't 

13406 

13447 

"r 

7/2 

13519 

13507 

1 3539 

“F 

9/2 

1L808 

U775  ' 

1 

14721 

15997 

15903 

15947 

“•G  , 

5/2 

17121 

17065 

17102 

2G 

7/2 

17288 

17196 

17228 

'•G 

7/2 

19057 

18959 

18910 

rms  deviation  (cm’M 

85.14 

33.34 

‘*Caicuiatod  by  using  parameters  of  W.  T.  Carnaii/  P.  R.  Fields, 
and  K.  Rjjnak  (CFR)  {J.  Chetn.  PhySs  , £9  (1968),  4412-65] , tabic  VI, 
line  2,  and  centered  to  experimental  centroids . 

^Calculated  by  using  derivatives  of  centroids  with  respect  to 
CFR  parameters  to  fit  CaWOu  centroids.  Exact  centroids  by  using 
parameters  table  VI,  lino  3,  differ  from  these  by  an  rms  of  0,27  , 

^Donald  £'.  Wortman,  Clyde  A.  Morrison , and  Nic)t  Karayianis , Harry 
Diamond  Laboratories  TR“1794  (June  1977). 


stark  splittings  for  Nd3'^:CaW04  by  using  the  crystal  field  parameters 
determined  by  the  original  fitting  procedure.  Deviations  between  the 
theoretical  energy  levels  of  the  two  calculations  were  1.8  cm'^  or  less 
for  the  lowest  53  Stark  levels  and  1.94  rms  for  the  highest  11  levels 
shown  in  table  VIII.  This  exercise  indicates,  fortunately,  that  the 
aqueous  eigenvectors  and  reduced  matrix  elements  are  adequate  for 
Nd^'^’crystal  field  calculations  at  least  to  the  degree  of  accuracy 
obtainable  by  the  seven-parameter  free-ion  Hamiltonian  considered  here. 
Because  of  the  inability  of  this  theory,  however,  to  obtain  a better  fit 
to  free-ion  centroids,  one  cannot  at  this  stage  predict  the  effects  on 
crystal  field  calculations  of  a more  profound  free-ion  theory. 


TABLE  VIII.  EFFECT  OF  Nd3+  FREE-ION  WAVE  FUNCTIONS  ON  THEORETICAL 
STARK  LEVELS  FOR  THREE  EXCITED  MULTIPLETS 


Multi plet 

Stark  energies  for  CaWO^  ( 

cm'l ) 

2u 

Theoret 1 ca 1 ^ 

Theoretical* 

lAl 

“G5/2 

3 

16979.7 

16980.9 

1.2 

(17101.8) 

1 

17041  .3 

17041.9 

0.6 

3 

17098.0 

17100.6 

2.6 

^<■7/2  (') 

1 

17223.1 

17225.5 

1 

2.4 

(17227.9) 

3 

17251.5 

17249.7 

1.8 

1 

17303.2 

17300.2 

3.0 

3 

17't26.6 

17423.1 

3.5 

’*<•7/2 

I 

18809.0 

18808.2 

0.8 

(18911.1) 

3 

18903.3 

18904.5 

1.2 

3 

18993.7 

18993.9 

0.2 

1 

19013.3 

19014. 1 

0.8 

rms  1 . 9^ 

3 

Free-ion  wave  functions  by  using  aqueous  psramotors 
of  W.  T.  Carnall,  P.  fi.  Fields,  and  K.  Rajnik  (CFR)  (J. 

Cliem.  Phys.f  £9  (1968),  4412-55],  table  VI,  line  2,  and 

CaWOi^  crystal  field  parameters  * 496»6b2;  » ] 

-867,308;  64,.  = 1032.  75;  B(.q  -8.09;  and  » 900.433  ' 

+ 1255.00  cm~^ . 

^Free-ion  wave  tunctions  by  using  fkiramj'tt'rs 

(table  VI,  line  3)  and  crystal  field  parameters  given  above.  ' 
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In  connection  with  the  best-fit  free-ion  parameters  for 
Nd^'*';CaW04  (table  VI,  line  3),  they  are  the  best-fit  parameters  only  to 
within  certain  limitations  imposed  on  their  deviation  from  the  aqueous 
parameters.  The  rigorous,  best-fit  parameters  gave  a fit  to 
28.23  cm"l  (instead  of  the  33.34  cm'^  fit  obtained  with  the  above 
parameters) , but  deviated  wildly  from  the  aqueous  values,  for 
example,  a,  and  y = -40.47,  727.92,  and  -13,284.7,  respectively. 
These  deviations  were  pared  to  deviations  compatible  with  the  parameter 
fluctuations  that  CFR**  found  over  the  lanthanide  series  and  then  were 
held  constant  while  best  fits  with  the  and  ^ were  obtained.  The 
result  of  this  exercise,  though  somewhat  arbitrary,  was  nevertheless 
physically  palatable  and  was  sufficient  for  performing  the  tests  on  the 
theory  described  above. 

Table  IX  shows  the  improved  free-ion  centroid  fit  obtained  for  the 
lowest  10  multiplets  of  Er^'^'iCaWO^  by  varying  the  free-ion  parameters. 
Similar  restrictions  on  a,  and  y to  those  described  above  were 
imposed  and  resulted  in  the  CaWO^  parameters  given  in  table  VI,  line  6. 
Perhaps  a more  liberal  leeway  on  the  a,  B#  and  y deviations  should  have 
been  allowed  in  this  fit  because,  whereas  the  rigorous,  best  fit 
parameters  gave  a fit  to  3.32  cm'^,  the  constrained  parameters  gave  a 
fit  only  to  47.73  cm"^.  Nevertheless,  the  CaWO^  free-ion  parameters 
obtained  with  the  constraint  differ  sufficiently  from  the  aqueous 
parameters  of  CFR**  to  perform  a test  of  the  theory  similar  to  that  for 
Nd  described  above. 

In  the  48  Stark-split  levels  in  the  10  lowest  multiplets  of  Er^'*'  in 
CaWOi,  calculated  by  the  two  sets  of  free-ion  wave  functions,  theoretical 
differences  were  only  1.13  cm"^  rms,  even  less  than  differences  found  in 
Nd^"*".  This  calculation  gives  further  strong  evidence  of  the 
insensitivity  of  crystal  field  calculations  to  free-ion  parameter 
changes  between  host  materials. 

As  a result  of  these  exercises,  we  conclude  that,  within  the 
limitations  of  the  seven-parameter  free-ion  Hamiltonian,  the  reduced 
matrix  elements  for  the  C in  the  crystal  field  Hamiltonian 


H = S B,  l.C,  (i) 
X km  km  1 km 


T.  Carnall , P.  R.  Fields,  and  K.  Rajnak,  J.  Chem.  Phys.,  49 
(1968),  4412-55.  — 
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TABLE  IX.  Er^-*-  MULTIPLE!  CENTROIDS  (crn'l) 


Mu  1 1 i p let 

1 

Carnal  1 
1 a 

1 parameters 

Best  fit^ 
parameters 

Experimental^ 
(CaWO^  centroids) 

15/2 

121. 

155 

13/2 

66<«5 

6681 

6632 

"I  11/2 

10255 

10258 

1021.1 

‘*'9/2 

121116 

121.09 

121.89 

“F9/2 

15280 

15321. 

15297 

‘**3/2 

181.97 

I8I.5I. 

181.03 

1/2(2) 

19293 

19190 

1911.6 

"F7/2 

20^457 

1 

20537 

20519 

“'^5/2 

22109 

22092 

22178 

"P3/2 

221.57 

221.89 

221.92 

rtns  deviation  (cm"l) 

68.10 

1.7.73 

Calculated  by  using  parameters  of  W.  T.  Carnall,  P.  R.  Fields, 
and  K,  Rajnak  (CFR)  [J.  Chem.  Phys.,  ^ (1968),  4412-55] , table  VI, 
line  5,  and  centered  to  experimental  centroids. 

^Calculated  by  using  derivatives  of  centroids  with  respect  to 
CFR  parameters  to  fit  Cah'Oi^  centroids . Exact  centroids  by  using 
parameters  of  table  VI,  line  6,  differ  from  these  by  an  rms  of 
3.42  cm~^ . 

^Donald  E.  Wortman,  Clyde  A.  Morrison,  and  Nick  Karayianis , Harry 
Diamond  Laboratories  TR-1794  (June  1977). 


1 


evaluated  in  the  aqueous  free-ion  basis  wave  functions  may  be  used 
independently  of  host.  These  reduced  matrix  elements  are  given  for 
Pr^'*',  Tm^'*",  Nd^'*' , and  Er^'*'  for  k = 2,  4,  6 in  tables  X to  XXI,  These 
tcibles  give  ( J' (L' S 'to ' ) | | E (i)  | | J (LSu) ) , * where  a general  crystal 
field  matrix  element  is  evaluated  according  to 


<J'M*  (L'S'u')  |E^Cj^(i)  |jM(LSo))>  = <J  (M)  k (m)  | J' (M' ) > 

>«  [ (L'S’to')  I |EXj^(i)  I I J(LSto))  . (15) 


*Multlply  these  ^reduced  matrix  elements , respectively , by  -*'15/28 , 
'/22/2S,  and  -*729/700— values  of  (<3(0)k(0)  \ 3(0)>*7)~'^—for  k = 2,  4, 

and  6 to  obtain  the  reduced  matrix  elements  (| \ |). 
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TABLE  X.  Pr3+  REDUCED  MATRIX  ELEMENTS  ( J ' ( L ' S ' ) | j EC2 ] ] J (LS)  ) 


n ' 

fl  ' 

n 

0 

” “ - 

(• 

n 

, 6 1.  ) 

n 

n 

n + 1 
0 

n *2 

(j 

n 

0 

n 

0 

n +5 
u 

n +6 
0 

1 

}P  0 

0.596362 

-0. 742609 

-0. 165856 

2 

IS  0 

-0.352898 

-0.085191 

-0.966218 

3 

3P  1 

3 

-0.51*7722 

-0.887462 

-0. 708172 

0.379780 

- 1 .032794 

k 

3P  2 

4 

0.855283 

-0.244479 

0.045818 

-0.693723 

-0.878429 

0.009562 

-1.121449 

3F  2 

0.339609 

0. 159808 

-0. 197953 

0.052284 

-0.974662 

-0.008547 

6 

ID  2 

6 

0.836991* 

0.2391*56 

1 .061070 

-0.070405 

'0. 7^^106 

7 

3F  3 

7 

0. 31*1565 

-0.216014 

-0.31*9532 

-0.088500 

-1 .083203 

8 

3F 

8 

0. 158500 

-0.186944 

0.385431 

'0.233227 

-I.02K30 

0. 4 14422 

9 

3H  it 

9 

-1.205814 

0.047739 

0.452078 

-0.014099 

-0. 131521 

10 

IG  L 

10 

-0.028045 

-0.2681 38 

-0.705772 

-0.666231* 

1 1 

3H  5 

1 1 

-1 . 309959 

0.449048 

0.024027 

12 

3H  6 

12 

-1 .520235 

0.080881 

13 

M 6 

13 

-3.027486 



L.  . 

. . ..  . 

The  (LSto)  are  the  assigned  Russell-Saunders  designations  for  the 
free-ion  wave  functions  according  to  the  discussion  in  section  3.  To 
save  space  and  avoid  printing  reduced  matrix  elements  that  are 
identically  zero  because  the  triangle  condition  among  J,  k,  and  J'  is 
not  met,  column  3 in  each  of  these  tables  gives  the  number  of  the  first 
multiplet  (n  ) to  which  the  multiplet  in  column  2 connects.  Across  a 
row,  the  reduced  matrix  elements  correspond  to  | 1 | as  n 
increases  in  steps  of  unity  from  the  value  rt  listed  in  column  3.  The 
lower  half  of  the  matrix  of  reduced  matrix  elements  is  given  by 

(j|  I J')  = (j- 1 IeCj^I  I j)  . (16) 


with  these  tables,  one  does  not  need  to  know  the  specific  structure 
of  cin  individual  free-ion  wave  function.  One  may  (1)  take  all  or  any 
particular  subspace  of  the  multiplets;  (2)  set  up  crystal  subspaces  of 
states  with  given  M that  are  connected  by  l!x;  (3)  evaluate  the 
interaction  matrix  according  to  equation  (15)  and  the  centroid  energy 
parameters  E(J(LSw)),  where 


TABLE  XII.  Pr3+  REDUCED  MATRIX  ELEMENTS  ( J ' (L ' S ' ) | | IXe | j J (LS) ) 


n ' 

n ' 

■ 

n 

o 

n = ' 

( 

, C.  t 

n n 

/"n 

o 

n + 1 
0 

fi  *2 
CJ 

.1. 

n ♦S  n *6 

0^0 

1 

3P  0 

12 

0. 3^i432S 

-0.0767bS 

2 

IS  0 

12 

0.01 7/89 

'0.9bl 341 

3P  1 

1 1 

-0. 381690 

'0.4SOS7b 

-0.029131 

It 

3P  ? 

8 

0. 129830 

0.47019s 

-0. 1661 39 

0.969251 

0.296966 

-C .982666 

6 

3F  2 

0. 382b06 

0.938258 

0.106357 

- 1 .037022 

0.703888 

-0.251998 

6 

ID  2 

8 

0. 186369 

-0.291392 

-0.361698 

0.028572 

-0. 103263 

- 1 ,6  72322 

7 

3f  3 

7 

-0. 31938S 

-0.092963 

1 .067999 

-0.291982 

6.000000 

-1.179963  -0.062869 

8 

3F  3 

8 

-0.268766 

0.889623 

0 . 76  36  7 1 

991 322 

0.862787 

0.896  305 

3 

Ji3  It 

9 

0.658622 

0.208596 

-0.998929 

0.977071 

-0. 197532 

10 

IG  It 

10 

-1.003319 

0.83668«t 

0.636236 

-1 .037987 

1 1 

3H  S 

1 1 

0.8Ii5I.’3 

-1.023579 

-0. 01.4  769 

12 

3H  6 

12 

1 . 1 39016 

0.086097 

13 

1 1 6 

13 

-0.950876 

TABLE  XIII.  Tm3+  REDUCED  MATRIX  ELEMENTS  ( J ' (L ' S ' ) j [ TC • | 1 J (LS) ) 


I 

! 

1 

n ' 1 

j 

1 

n = ^ 

(t 

^ J 

ft  + 1 n *? 

«.  ' 0 

I r,  .3 

n 1 9 1 

“ i 

1 n 

\ - - *'* 

11  +6 
0 

1 , 

3P  0 1 
1 

-0.978731  i 

0.812198  J 

-0.229881  1 

1 

\ 

7 

1 

;s  0 : 

1 

9 

-0.766095 

-0.305996 

0.629198 

1 

3 1 

, 3P  1 1 

3 

1 

0.597722  1 

0.679662  ! 

0.512369 

0.916868 

1 .032799 

1 

4 i 

3P  2 ' 

9 

1 

-0.8169  50  ^ 

0.089385 

0.026101 

0.513109 

0.993808 

0. 71 3875 

1 -1.069602 

s 

3f  2 ' 

6 I 

-0.510792  i 

0.395813 

0.085387 

-0.797708 

0.799068 

0. 109093  1 

1 

£> 

10  2 ^ 

6 i 

-0.605753 

0.559295 

1.036016 

0.987717 

0.586719  ' 

; 

3P  3 

7 

-0.391666 

0.067899 

0.591389 

, 0.137911 

1 .083203 

8 

3F  4 1 

i ^ 

-0.155608 

-0.991290 

* 0.081259 

0.912899  1 

1 1.001638 

0.  337688 

9 

' 3H  4 

9 

0.706729 

0.590192 

-0.156619 

! 0.665526 

1 

, -0. 599806 

to 

IG  4 

10  ' 

0.609239 

0. 568366 

-0. 300152 

O.6330IS 

J 1 

! 3H  9 

1 1 

1 . 509959 

-0.997752 

1 

, 0. 04)^24 

12 

1 5H  6 

12  ' 

1.529082 

0.190710 

13 

1 1 6 

13 

1 3 . 0 1 86  58 

i 

1 

1 

TABLE  XV.  Tm3+  REDUCED  MATRIX  ELEMENTS  ( J ' (L  ‘ S ' ) | j r,C(,  j ] J (LS)  ) 


M 
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■ *2 

,,  *3 

n *4 
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n 

fj  *6 

0 

12 
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.1-  1 72/00 

• 

1 

O.OIS^21 

'■.■*^6233 

i 

I 

t 1 

. 33I69C 

--8:1688 

'0.0‘«?'03 

m 2 

8 

- 1 MOISI... 

'0.033Ib2 

■0.316036 

-0.853503 

-0. I9062j 

-1.07:767 

2 

8 

-0  38'i6  33 
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0.979178 

-0.688925 

-0. 783873 

b 

ID  2 

8 

■ ■ - l^).'0  7l 

-.7.603636 

-0.037326 

0. 1631 73 

-0.3891:9 

1 . 169055 

/ 

if  3 

7 

0. 

■'.87/810 

-0.681971 

-0.699555 

-0.000000 

:. 171510 

-0. 109696 

<i 

3‘  - 

8 

0.662  1 1 I 

-0.679100 

0-383268 

-1.230376 

-0.623888 

0.798383 

a 

I* 

0.  (788  70 

-(>■  776827 

0. 123802 

-0.988090 

-0.398150 

10 

IG  ^ 

10 

0.4:632'- 

-..9388  76 

0.  182  788 

1 .02001 1 

1 1 

3-^  i 

1 1 

-j.—iiii 

1 .020878 

-0.096663 

i- 

5h  b 

12 

- 1 . i2t*;i2 

0. 187989 

n 6 

13 

0.h4^',,’2 

. ... 

and  (4)  vary  the  B)^j^  and  centroids,  E(J(LSto)),  to  obtain  a least  squares 
fit  with  experimental  levels.  Within  the  approximation  of  constant 
basis  wave  functions,  this  procedure  is  tantamount  to  a simultaneous 
diagonalization  of  H + H since,  by  Wong's  method,  the  "free-ion" 
Hamiltonian  may  be  inferred  from  the  resultant  centroids,  and  the 
crystal  field  Hamiltonian  is  determined  by  the  resultant  B 

km 
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TABLE  XIX,  Er3+  REDUCED  MATRIX  ELEMENTS  { J ' (L ' S ' oj' ) 1 | ECj | | J (LSu) ) 


n' 

n 

0 
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1 UC2I 

n 

0 

n +1 
0 

n +2 
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9 
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10 
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10 
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1 
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TABLE  XX.  Er3+  REDUCED  MATRIX  ELEMENTS  (J  ' (L  ’ S ' w ' ) | | ZC4  | ] J (LSu)) ) 
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TABLE  XXI.  Er3+  REDUCED  MATRIX  ELEMENTS  { J ' ( L ' S ' oi ' ) | j ECg | 1 J ( LSw) ) 
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